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ABSTRACT 

We present the results from the spectral analysis of more than 7,500 RX TE spectra 
of 10 AGN. Our main goal was to study their long term X-ray spectral variability. 
The sources in the sample are nearby, X-ray bright, and they have been observed by 
RXTE regularly over a long period of time (~7-ll years). High frequency breaks have 
been detected in their power spectra, and these characteristic frequencies imply time 
scales of the order of a few days or weeks. Consequently, the RXTE observations we 
have used most probably sample most of the flux and spectral variations that these 
objects exhibit. Thus, the RXTE data are ideal for our purpose. Fits to the individual 
spectra were performed in the 3-20 keV energy band. We modeled the data in a 
uniform way using simple phenomenological models (a power-law with the addition 
of Gaussian line and/or edge to model the iron Ka emission/absorption features, if 
needed) to consistently parametrize the shape of the observed X-ray continuum of the 
sources in the sample. We found that the average spectral slope does not correlate with 
source luminosity or black hole mass, while it correlates positively with the average 
accretion rate. We have also determined the (positive) "spectral slope - flux" relation 
for each object, over a flux range larger than before. We found that this correlation is 
similar in all objects, except for NGC 5548 which displays limited spectral variations 
for its flux variability. We discuss this global "spectral slope - flux" trend in the 
light of current models for spectral variability. We consider (i) intrinsic variability, 
expected e.g. from Comptonization processes, (ii) variability caused by absorption 
of X-rays by a single absorber whose ionization parameter varies proportionally to 
the continuum flux variations, (iii) variability resulting from the superposition of a 
constant reflection component and an intrinsic power-law which is variable in flux but 
constant in shape, and, (iv) variability resulting from the superposition of a constant 
reflection component and an intrinsic power-law which is variable both in flux and 
shape. Our final conclusion is that scenario (iv) provides the best fit to the data of all 
objects, except for NGC 5548. 
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1 INTRODUCTION 

The current paradigm for the central source in AGN pos- 
tulates a black hole (BH) with a mass of 10 6 — 10 9 Mq, 
and a geometrically thin, optically thick accretion disc that 
presumably extends to the innermost stable circular orbit 
around the black hole. This disc is supposed to be respon- 
sible for the broad, quasi-thermal emission component in 
the optical-UV spectrum of AGN (the so-called Big Blue 
Bump). At energies above ~2 keV, a power-law like com- 
ponent is observed. This is attributed to emission by a hot 
corona (T ~ 10 s - 10 9 K) overlying the thin disc. The corona 
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up-Comptonizes the disc soft photons to produce the hard 
(E ~ 2 — 100 keV) X-ray emission. An Fe Ka line is also 
detected, produced either by fluorescence on cold matter at 
E — 6.4 keV or by photoionization of H-like or He-like Fe at 
6.9 keV and 6.7 keV, respectively. This feature is believed to 
be emitted on the surface of the accretion disc, from the re- 
processing of the X-rays. In many cases the line is resolved, 
and its width can, in theory, constrain the position of the 
inner radius of the disc and the spin of the central BH. 

The AGN X-ray emission is strongly variable, even on 
time scales as short as a few hundred seconds, both in flux 
and in the shape of the observed energy spectrum. These 
variations can provide information on the physical condi- 
tions, the size, and the geometry of the X-ray emitting re- 
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gion. Regarding the flux variability, although the broad band 
power spectra of AGN are not well understood, the identifi- 
cation of characteristic time scales in them, combined with 
the large range in the masses of the BHs that have been stud- 
ied so far, can constrain (to some d egree) the size and th e 
geometry of the X-ray source (e.g. lM c Hardv et al.l l2006). 
On the other hand, recent detailed studies of the spectral 
variations observed in a few objects provide valuable de- 
tails of the physical conditions of the X-ray reprocessing 
|Miniutti et alj|2007t see e.g.) and/or on the physical prop- 
erties of absorbing m aterial which may exist close to the 
central source (see e.g. iMiller. Turner fc Reeves! [2008 ), 

Progress in settling alternative interpretations of the 
observed spectral variability can be made by correlating the 
spectral and t iming properties in AG N. Along this line, in a 
recent paper (|Papadakis et al]|2009l ) we studied 11 nearby, 
X-ray bright Seyferts, which have been intensively observed 
over the last two decades by the Rossi X-ray Timing Ex- 
plorer (RXTE). We detected a positive correlation between 
their mean spectral slope and the characteristic break fre- 
quency, Ubi, at which the power spectrum changes its slope 
from -1 to -2. To the extent that this frequency corresponds 
to a (still unknown) characteristic time scale of the X-ray 
source, this correlation suggests that the mean X-ray spec- 
tral slope in AGN is determined by a fundamental parameter 
of the system, most probably its accretion rate. If true, this 
result has important implications on current Comptoniza- 
tion models. 

In this work we use the results from the model fitting 
of thousands of RXTE spectra of 10 AGN t hat were in- 
cluded in the sample of iPapadakis et ail l|2009l) . Our prime 
goal is to study in detail the X-ray spectral variability of 
each object. Although the spectral resolution of RXTE is 
not as good as that of Chandra, XMM-Newton or Suzaku, 
the RXTE data sets that we used are ideal for our purpose. 
First, the sources in the sample are nearby and X-ray bright. 
As a result, their X-ray spectral shape can be accurately de- 
termined even from relatively short (~1— 2 ksec) RXTE ob- 
servations. Second, each source in the sample was observed 
regularly (more than 200-300 times, and in a few cases more 
than a thousand times) over a long period of time (of the 
order of 7-11 years). This is contrary to typical Chandra 
and XMM-Newton observations of individual AGN, which 
usually span a period of a few days (maximum). Third, the 
characteristic frequencies of the objects imply time scales of 
the order of a few days or weeks. Consequently, the RXTE 
observations we have used most probably sample most of 
the flux and spectral variations that they exhibit. 

A subsamp le of the observat i ons w e study have already 
been studied bv lPapadakis et alj (|2002l ) . These authors used 
hardness ratios in various energy bands to study the spectral 
variability of four objects, namely NGC 5548, NGC 5506, 
NGC 4051 and MCG -6-30-15. In the present work we stud- 
ied many more observations of these objects and we have 
also analysed data for six more Seyferts. Furthermore, in- 
stead of computing hardness ratios, we fitted each spectrum 
with the same phenomenological model of a "power law plus 
a Gaussian line" to derive the observed photon index, r o b s - 
In this way, we parametrized all available 3-20 keV RXTE 
spectra of the sources in a uniform way. We then used the 
resulting best model fitting values to: i) determine the mean 
spectral slope for each object and investigate its correlation 



with the main source parameters such as BH mass, accre- 
tion rate and luminosity, and ii) construct their long term 
"spectral slope vs. flux" relations. For the reasons explained 
above, we believe that we determined the flux related spec- 
tral variability in these sources more accurately than before. 

Several possible reasons for the observed X-ray spec- 
tral variability in AGN have been proposed in the 
past. The simplest possibility is that the r o b s vari- 
ations co rrespond to intrinsic variations in continuum 
slope (e.glHaardt. Maraschi fc Ghisellinilll997l;ICoppilll999l; 



iBeloborodovl Il999l ). For example, Nandra fc Papadakis! 
|200ll ) have showed that the X-ray spectral changes in 
NGC 7469 are intrinsic and respond to the UV-soft pho- 
ton input variations, exactly as one would expect in the 
case of thermal Comptonization models. However, there 
has also been certain observational evidence that the F b s 
variations originate as a result of superposition of differ- 
ent spectral components with different variability proper- 
ties. One such possibility is a spectrum composed of a 
constant reflection component and a variable, in normal- 
izatio n only, power-law like con t inuum of constant slope 
(e.g. iTavlor. Uttley fc M c Hardvl 120031 ; iPonti et al l 120061 ; 
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iMiniutti et al.l l2007. and references therein). Another possi- 



bility is a constant slope power law which varies in flux and 
complex variations in an absorber, e.g. its ionization state 
an d /or the covering facto r of the source (see recent review 
bv lTurner fc Miller! |2009| ). 

One of our main result is that the X-ray spectral slope 
steepens with increasing flux in a similar way for all the 
objects in the sample. This similarity implies that, although 
all the above mentioned mechanisms may operate in AGN, 
it is just one of them which is mainly responsible for the 
observed spectral variability in AGN. 

The paper is organized as follows. In Section 2, we dis- 
cuss the data selection and reduction, and we also describe 
the phenomenological models that we used to parametrize 
the shape of AGN continuum. In Section 3, we report the 
spectral fitting results, and the results from the spectral vari- 
ability study. In Section 4, we discuss our findings and com- 
pare them with predictions of several spectral variability 
models. We give our conclusions in Section 5. 



2 DATA REDUCTION AND ANALYSIS 

2.1 Data selection and reduction 

In Tab. [T] we list the 10 AGN in our sample, some 
of their properties (distance, redshift and BH mass esti- 
mates), and details of th e RXTE observations we used in 
this work. We used the iPeterson etaf] l|2004l) estimates, 
from reverbaration mapping, for the BH mass (Mbh) in 
FairaU 9, NGC 4051, NGC 3227, NGC 5548, NGC 3783, 
and NGC 3516. In the case of Mrk 766 and MCG-6-30- 
15 we used th e stell ar vel ocity dispecrsion meas urements 
of iBotte et ail (|2004|) and iM^Hardy et all (|2005l ). respec- 
tively, and the iTremaine et al.l "( 20021 ) Mrh — a relation to 
estimate BH mass. The BH mass estimates for Ark 564 
and NGC 5506 were t aken from lZhang fc Wand (|2006l ) and 
IWang fc~Zh ang (200?]), resp ectively. They were also based 
on the Tremaine et al.l (|2002l ) Mbh — a relation, although in 
this case the stellar velocity dispecrsion was obtained from 
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Table 1. Summary of the RXTE observations and properties of the AGN in the sample. In column 2, N is the total number of 
observations considered for a given AGN. Dates of the first and last observation are listed in column 3. The luminosity distances, D, and 
redshifts, z, listed in columns 4 and 5, were taken form the NED database. In column 6 we list black hole mass estimates. 



Source 


N 


Start/End 


Lum. Distance 


Redshift 










D 1 Mpc 


z 


xlO 6 Mq 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Fairall 9 


671 


1996-11-03/2003-03-01 


199 


0.047 


255 


Akn 564 


505 


1996-12-23/2003-03-04 


98.6 


0.025 


1.9 


Mrk 766 


219 


1997-03-05/2006-11-07 


57.7 


0.013 


3.6 


NGC 4051 


1257 


1996-04-23/2006-10-01 


12.7 


0.002 


1.9 


NGC 3227 


1021 


1996-11-19/2005-12-04 


20.4 


0.004 


42.2 


NGC 5548 


866 


1996-05-05/2006-11-06 


74.5 


0.017 


67.1 


NGC 3783 


874 


1996-01-31/2006-11-08 


44.7 


0.010 


29.8 


NGC 5506 


627 


1996-03-17/2004-08-08 


29.1 


0.006 


28.8 


MCG-6-30-15 


1214 


1996-03-17/2006-12-24 


35.8 


0.008 


6.8 


NGC 3516 


250 


1997-03-16/2006-10-13 


37.5 


0.009 


42.7 



Table 2. Quality of the various model fits to the data. The number of model-fit "successes", n s , and global propabilities, P c (see Sec. 
2.2 for details), are listed for various models. 



PL PLG ePLG™ ePLG 6 ePLG c 

Source n s P c n a P c n a P c n a P c n s P c 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 



Fairall 9 


642 


0.81 










Akn 564 


486 


0.92 










Mrk 766 


208 


0.54 










NGC 4051 


1167 


5 x 10~ 4 


1199 


0.75 






NGC 3227 


919 


< 10~ 5 


980 


0.94 






NGC 5548 


802 


1 x 10~ 3 


827 


0.77 






NGC 3783 


635 


< 10~ 5 


799 


< io- 5 


821 


0.09 


NGC 5506 


382 


< 10~ 5 


556 d 


< io- 5 < d 


586 d 


0.05 d 


MCG-6-30-15 


925 


< 10- 5 


1076 


< 10~ 5 


1136 


0.02 


NGC 3516 


160 


< io- 5 


221 


< 10" 5 


227 


0.004 



1147 0.22 

227 0.004 231 0.05 



a the edge energy fixed at 7.1 keV 

b the edge energy allowed to vary 

c both the line and edge energies allowed to vary 

d PLG or ePLG model plus cold absorption with Nu = 3.7 X 10 22 cm~ 2 



the FW HM of the Qu i line, using the FWHM Om - a re- 
lation of[G rccne h Hoi (|2005l ). Since the reverbaration map- 
ping estimates have been obtained by "forcing" the AGN 
Mbh — <y relationship to have the same norm alization as the 
Mbh — c relationship for quiescent galaxies (|Peterson et al.l 
|2004| ). we belive that the usege of two different Mbh deter- 
mination methods should not add significant scatter to the 
correlations we present below. 

For each object we considered all the RXTE observa- 
tions that were performed until the end of 2006. Some of 
these objects were monitored with RXTE for as long as 
10 years. The reason we chose to study these AGN is pre- 
cisely because they have been observed frequently by RXTE. 
They are all nearby, X-ray bright, Type I objects (except for 
NGC 5506, which is optically classified as Type II Seyfert). 
Their central BH mass ranges from ~2xl0 6 M to ~2xl0 8 
Mq, and their accretion rate from a few percent of the Ed- 
dington limit (in NGC 3227) to almost close to this limit 
(in Ark 564; see e.g. M c Hardy et al. 2006). Consequently, 
they could be considered as representative of the Seyfert 
population in the local Universe. 

We used data from the Proportional Counter Array 



(PCA; Ijahoda et al.l Il996j ). The typical duration of each 
observation was ~l-2 ksec. The data were reduced using 
FT0QLS v.6.3. The PCA data were screened according to the 
following criteria: the satellite was out of the South Atlantic 
Anomaly (SAA) for at least 30 min, the Earth elevation an- 
gle was ^ 10°, the offset from the nominal optical position 
was < 0°02, and the parameter ELECTRON-2 was sC 0.1. 

We extracted STANDARD-2 mode, 3-20 keV, Layer 1, 
energy spectra from PCU2 only. For background subtraction 
we used the appropriate background model for the faint ob- 
ject^]] available from RXTE data analysis web-pages. The 
background subtracted spectra were rebinned using grppha 
so that each bin contained more than 15 photons for the 
X 2 statistic to be valid. We considered only these data sets 
for which binning resulted in at least 15 PHA channels. We 
used PCA response matrices and effective area curves cre- 
ated specifically for the individual observations by pcarsp. 



pca_bkgd_cmf aint 17_eMv2005 1128. mdl 
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2.2 Modeling the spectral shape of the AGN 

Fits to the individual spectra were performed in the energy 
range 3-20 keV, where the PCA is most sensiti ve, using the 
XSPEC v. 11.3.2 software package (|Arnaudlll996h . Our prime 
aim was to parametrize in a uniform way (i.e. by fiting the 
same model to all the observations) the shape of the ob- 
served X-ray continuum of the sources in our sample. We 
first fitted a simple power-law model (hereafter PL) to all 
the individual spectra of all sources. In the cases when the 
global goodness of the PL fit was unacceptable, we added 
additional spectral components. The form of these compo- 
nents was chosen based on visual examination of the PL fit 
residuals in the cases when the model did not fit the data 
well. 

We did not consider the Galactic absorption effects as 
the Galactic column to all sources is so low that cannot 
affect the spectrum above 3 keV. In the case of NGC 5506, 
all models were modified by cold absorption using wabs in 
XSPEC, with the hydrogen column density fixed at Nh = 
3.7 x 10 22 cm -2 l|Perola et al.ll2002l ). 

We assumed that a model fits well the spectrum of an 
individual observation if the probability of accepting the null 
hypothesis, p nu n, is larger than 0.05. Given that the number 
of observations for each AGN is quite large, it is possible 
that Pnun < 0.05 in some cases, even if the model is the 
correct one. Therefore, the global goodness of a model fit to 
all the spectra of each source had to be judged in a statistical 
way. One method to achieve this is to consider the fit to the 
individual spectra as a "success" when p nu n > 0.05, and as a 
"failure" otherwise. Under the hypothesis that the model is 
the correct one, the probability of a "success" ( "failure" ) is 
p s = 0.95 (l-p s = 0.05), by definition. In this approach, the 
model fit to each spectrum can be considered as a Bernoulli 
trial. 

For a given source and a given model we recorded the 
total number of spectra, N, the number of model-fit "suc- 
cesses" , n s , and the number of model- fit "failures" , rif . Then, 
we used the binomial distribution to estimate the probabil- 
ity, P c , that the number of "failures" will be larger than m in 
N trials by chance. This is equal to the probability that the 
numbers of "successes" will be smaller than n a . Therefore, 
P c can be estimated as 

p c = E n)d(i-Psf- j . a) 

J=0 V 3 / 

We accepted that a model fits globally the spectra of a source 
if 0.05 < P c < 0.95. We report below the results from the 
various model fits to the data of all sources. 



2.2.1 "Power-law" model fits 

The PL model fited well the spectra of only three sources: 
Fairall 9, Akn 564 and Mkn 766 (n B and P c values are listed 
in the second and third columns of Tab. [2). In the remaining 
7 objects we clearly detected 6-7 keV residua indicative of 
an iron K a line (see Fig. [TJ. 



2.2.2 "Power-law plus line" model fits 

Motivated by shape of the PL best model fit residua around 
6-7 keV, we added a Gaussian component to account for the 
Ka iron line in these spectra (powerlaw+gaussian in XSPEC, 
hereafter PLG). In all cases, the iron line energy and width 
were fixed at 6.4 keV and 0.1 keV, respectively. The goodness 
of the new model fits improved and became acceptable in 
NGC 4051, NGC 3227 and NGC 5548 (n s and P c values for 
the PLG model are listed in columns 4 and 5 of Tab.[2J). 



2.2.3 "Power-law plus line plus edge" model fits 

Figure [2] shows the best PLG model fit residua in the case of 
an unacceptable model fit to a spectrum of (a) NGC 5506, 
(b) NGC 3783, (c) MCG -6-30-15, and (d-e) NGC 3516. 
The best-fit model residuals indicate the presence of an 
edge-like feature around ~7 keV. We therefore considered a 
edge* (powerlaw+gauss) model in XSPEC (ePLG hereafter). 
We kept the edge energy fixed at 7.1 keV (the threshold en- 
ergy for neutral iron) and we re-fitted the individual spectra 
of these sources. 

The global fit became acceptable in the case of 
NGC 5506 and NGC 3783. The optical depth of the edge, 
t, was detected at the 2a significance level in 298 out of 627 
observations, in the case of NGC 5506, and in 319 out of 874 
observations in NGC 3783. The respective n s and P c values 
listed in columns 6 and 7 of Tab. [2] 

In MCG-6-30-15 the model fit improved significantly, 
and the model became acceptable, when we allowed the en- 
ergy of the edge to vary between 6 and 10 keV (results are 
listed in columns 8 and 9 of Tab.[2]|. The edge was detected 
at the 2a significance level in 863 out of 1214 data sets. Us- 
ing the best fit results for these spectra only, we estimated 
that -Bodge ~ 7.4 keV, an energy that is indicative of absorp- 
tion by ionized iron. However, the standard deviation of the 
best fit edge energy values is substantial (a ~ 1 keV) , which 
suggests that one should treat this result with caution (see 
discussion in Section 2.2.5). 

In the case of NGC 3516 we did not record an improve- 
ment to the global fit of the ePLG model, even when the 
edge energy was allowed to vary (compare columns 6 and 7, 
and 8 and 9 of Tab. [2]). Thus, we thawed the energy of the 
iron line and allowed it to vary between 5.5 and 7.5 keV. 
This extra modification provided a statistically acceptable 
global fit (the n s and P c values listed in columns 10 and 11 
in Tab. [2] correspond to this model). The average line en- 
ergy, and its standard deviation, calculated for the data sets 
where the line is detected at the 2a significance level (117 
of 250) were 6.1 and 0.3 keV, respectively. The absorption 
edge was detected at the 2a significance level in 122 out 
of the 250 spectra of this source. Using the best fit results 
of these observations only, we calculated that the average 
energy of the absorption edge, and its standard deviation, 
were -B c d gc = 7.9 and a = 1.1 keV, respectively. 

2.2.4 Summary of the acceptable model fits 

1. A simple PL model describes well the spectra of Fairall 9, 
Akn 564 and Mkn 766. 

2. In NGC 4051, NGC 3227, and NGC 5548 a statistically 
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NGC4051 
p=0.0264 



MCG-6-30-15 
p=0.0000 



NGC 3227 
p=0.0372 



NGC 3783 
p=0.0239 



NGC 5548 
p=0.0033 



NGC 5506 
p=0.0032 



NGC 3516 
p=0.0171 



10 



Energy (keV) 



(a) 



(b) 



(c) 



(d) 



(e) 



(f) 



(g) 



20 



Figure 1. Contribution to x 2 for the PL model fits to the data of 
(a) NGC 4051 (obsid 70141-01-01-00), (b) MCG-6-30-15 (40155- 
01-02-02), (c) NGC 3227 (50153-07-25-20), (d) NGC3783 (40152- 

04- 65-00), (e) NGC 5548 (30220-01-05-00), (f) NGC 5506 (50153- 

05- 67-10), and (g) NGC 3516 (40223-01-84-00). The fit probabil- 
ities are indicated on each panel. Residuals indicative of an iron 
Ka line are visible in the 5-7 keV band. 



acceptable global fit to the data is obtained with the PLG 
model. 

3. In NGC 3783 and NGC 5506, the addition of an absorp- 
tion edge with energy fixed at 7.1 keV is needed to achieve 
an acceptable model fit (ePLG). 

4. In MCG 6-30-15 the best global fit is reached by the 
ePLG model but only when the edge energy is allowed to 
vary (the resulting mean edge energy is -E c d g o ~ 7.4 keV). 

5. Finally, the same ePLG model provides an acceptable 



NGC 5506 
p=0I)Q84 



(a) 



NGC 3783 
p=0.0015 



(b) 



MCG-6-30-15 
p=0.0148 _L 



(c) 



NGC 3516 
p=0.0357 



(d) 



_L 



1 •— 

NGC 3516 
p<10 4 



(e) 



10 



20 



Energy (keV) 



Figure 2. Contributions to x 2 for the PLG best fit models ap- 
plied to exemplary data of (a) NGC 5506 (obsid 90145-01-01-00), 
(b) NGC 3783 (obsid 92113-03-81-00), (c) MCG-6-30-15 (obsid 
20310-01-01-01) and (d-e) NGC 3516 (obsids 60131-09-02-00 and 
20316-01-02-00). The fit probabilities are indicated on each panel. 
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Figure 3. The average r o b s (a) and EW of the iron line (b) 
in NGC 3516 for different spectral models. The model labels are 
Ml - PL; M2 - PLG; M3 - ePLG; M4 - same as M3 but with 
variable edge energy, E e( ^ se ; M5 — same as M3 but with variable 
line energy, En nE ; M6 — same as M3 but with variable both E ca \ sa 
and 7?ii no . The error- bars represent standard errors on the mean. 
The horizontal lines indicate the average r o b s and EW. 
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global fit to the NGC 3516 spectra but only when both the 
energy of the iron line and the energy of the absorption edge 
were left free to vary (Eu nc ~ 6.1, and -E e d gc ~ 7.95). 

In Section 3, we report the results from: (a) the study 
of the correlation of the average spectral slope with other 
source parameters, such as luminosity, accretion rate and 
black hole mass, (b) the study of the spectral slope variations 
as a function of the source luminosity for each object in the 
sample, and (c) the study of the average line's equivalent 
width variations with other source parameters. In each case, 
we use the appropriate, for each source, best fitting slope, 
r o bs, and line's EW values. 

We do not consider the implications from the "variable 
edge energy" model fits to MCG -6-30-15 or from the "vari- 
able iron line and edge energy" model fits to NGC 3516. We 
believe that the Sedge and i?ii ne best fit values from these 
model fits should be treated with caution due to the low 
energy resolution of the PC A instrument. The residua plot 
in Figs. [2]i-e show this clearly in the case of NGC 3516. It 
is the line-like positive residual at ~5 keV that forces the 
best fit iron line energy to be less than 6.4 keV in this case. 
However, we can also get a good fit to the the spectrum from 
e.g. the 20316-01-02-00 observation of NGC 3516 (Fig. [2^) if 
we consider the ePLG model with -Eii ne and -E c d gc kept fixed 
at 6.4 keV and 7.1 keV, respectively, together with intrinsic 
neutral absorption, i.e. if we add a wabs component to the 
model and we let Nh to be a free model parameter. 

We note that in the case of NGC 5506, NGC 3783, 
MCG 6-30-15 and NGC 3516, the F obs values were not dra- 
matically affected when we added an absorption edge or 
when we let -E e d ge or -Eime to vary freely during the model 
fitting process. For example, Fig. [3^ shows the mean r o b a for 
NGC 3516 in the case of all 6 models that we applied to its 
individual spectra (PL, PLG, ePLG and its modifications: 
ePLG model with variable E a dg C and/or En nc ). Despite the 
fact that, statistically speaking, models 1-5 are not formally 
acceptable, the mean best fit spectral slope values are almost 
identical. The largest difference is the one between F b s in 
the case of the PL model (Ml in Fig. [3} and in the case of 
the ePLG model when both the edge and line energies are 
left free to vary (M6 in Fig. [3} . Although the difference be- 
tween the mean T's is statistically significant (the error-bars 
show standard errors on F), it is so small (~ 0.1) that it 
cannot affect any of the results we discuss in Section 3. 

We reached the same conclusion for the mean EW of the 
iron line. In Fig.[3p, we plot EW in the case of the 5 models 
that include an emission line component. The averages were 
calculated over the data-sets in which the line was detected 
at 2a level. All five measurements are in agreement with 
each other (again, the error-bars show standard errors on 
EW). The largest difference is the one between the average 
EW in the case of model M4 and model M6 (~100 eV). 
Taking into account their errors, the two measurements are 
consistent with each other, and so this difference does not 
affect any of the conclusions we present in Section 3. 



3 SPECTRAL VARIABILITY ANALYSIS 

Figure [4] shows r o b s and the X-ray luminosity, L2-10 plot- 
ted as a function of time for NGC 4051, MCG-6-30-15, 



NGC 3516 and Mrk 760 The first two objects are represen- 
tative of sources observed by RXTE frequently (on average 
every 2 — 3 days), between 1996 and 2006. On the other hand, 
NGC 3516 and Mkn 766 were systematically observed by 
RXTE only between the years 1999-2002 and 2004-2007, 
respectively. Outside of these dates data for these sources 
are scarce. Their light curves are representative of sources 
which were observed less than ~700 times by RXTE from 
the start of the mission until the end of 2006. 

The (unweighted) mean photon indices, F b s , and 2- 
10 keV luminosities, L2-10, for the sources in our sample 
are listed in Tab. [3] (the standard error on F b s and -L2-10 
are very small, mainly because of the large number of data 
points). Application of the standard \ 2 test to the r o b a and 
£2-10 light curves suggests that all sources show significant 
spectral and flux variations (reduced \ 2 values are listed 
in Tab. [3] the number of data points is large enough that 
even in the case when % 2 /d.o.f. ~ 1.3-1.4 the variations we 
detected are significant at more than the 3<r level). 

To further quantify the variability properties of the 
AGN in our sample, we calculated the excess variance, a 2 s , 
of the Fobs and L2-10 light curves as follows, 

^xs = y^(Xi-x) 2 yVerr.i- (2) 

i=l i=l 

The excess variance is a measure of the scatter about 
the mean in a light curve, corrected for the contribution ex- 
pected from the experimental errors. In the formula above, 
x stands either for r o b s or L2-10, and <J 2 rrii is the arith- 
metic mean of the upper and lower 90% confidence error on 
Fobs and Z/2-io- We then estimated the fractional root mean 
square variability amp litude, F va , T = \/ai B /x 2 and its error 
using equation (B2) in lVaughan et afT (120031 ') (this error ac- 
counts only for the uncertainty due to the experimental error 
in each point). 

The results are listed in Tab. [3] In terms of variabil- 
ity amplitude, the flux variations are stronger (i.e. of larger 
amplitude) than the spectral variations. Furthermore, there 
seems to be a tendency for the objects with stronger flux 
variations to show stronger spectral variations as well. To 
quantify this correlation, we used the Kendall's r test. We 
found that r = 0.51, which implies that the positive correla- 
tion between F var - and F var l( 2 _io) ^ s significant at the 4% 
level. 

3.1 Correlation between the average spectral 
slope and other source parameters 

Figure [5] shows the average observed photon index, r o b s , 
plotted as a function of (a) BH mass, (b) average 2-10 lu- 
minosity, and (c) average X-ray mass accretion rate in Ed- 
dington units, rax,E = Lz-io/Le, where Le = 1.3 x 10 38 
M/Mq ergs s -1 is the Eddington luminosity of an AGN 
with a black hole mass M. This quantity should be repre- 
sentative of the total accretion rate (i.e. Lbd/L^) of each 

2 The luminosity is estimated as L2-10 = 47r_D 2 F2_io, where 
D is the luminosity distance, listed in Tab. [T] and F2— 10 is the 
observed 2-10 keV flux, from the appropriate best model fitting 
results for each source. The flux measurements are corrected for 
absorption in the case of NGC 5506. 
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Table 3. Variability properties of r o b s and L2-10 based on the best-fit model parameter values. 
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Target 



Fobs 



X 2 /d.o.f. 



(%) 



£2-10 

(xlO 42 ) erg s" 



X 2 /d.o.f 



(%) 



Fairall 9 


1.794±0.006 


1.3 


3.7±0.3 


102±1 


16 


26.0±0.4 


Ark 564 


2.48 ±0.01 


1.6 


5.0±0.3 


25.5±0.4 


20 


30.7±0.3 


Mrk 766 


2.02 ±0.01 


2.3 


4.6±0.4 


10.6±0.2 


55 


30.3±0.5 


NGC 4051 


1.842±0.009 


6.9 


14.9±0.3 


0.493±0.007 


83 


49.7±0.2 


NGC 3227 


1.552±0.008 


7.1 


13.7±0.2 


1.86±0.03 


109 


48.0±0.2 


NGC 5548 


1.735±0.004 


2.1 


2.9±0.2 


28.8±0.4 


196 


39.4±0.2 


NGC 3783 


1.616±0.004 


2.5 


4.7±0.1 


17.0±0.1 


111 


23.2±0.1 


NGC 5506 


1.800±0.003 


2.4 


3.1±0.1 


10.2±0.1 


169 


22.79±0.09 


MCG-6-30-15 


1.731±0.005 


3.6 


7.2±0.2 


7.48±0.06 


93 


27.2±0.1 


NGC 3516 


1.40 ±0.01 


7.6 


14.2±0.5 


6.5±0.2 


179 


41.3±0.3 
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Figure 4. The RXTE light curves of r obs and L2-10 for objects with the best (NGC 4051 and MCG-6-30-15, a-b) and worst (NGC 
3516 and Mrk 766, c-d) time coverage. 
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Figure 5. Relations between the average photon index, r o t, s , and (a) BH mass, (b) average 2-10 keV luminosity, and (c) average mass 
accretion rate in Eddington units. The correlation in (c) is best described with a power-law, r o b s <x m x E ' (solid curve). 
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Figure 6. Relations between the average EW of the iron line and the average (a) 2—10 luminosity, (b) mass accretion rate, and (c) 
photon index. The solid line in panel (a) indicates the best fit relation EW oc L 2 _ io 
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Figure 7. Dependence of r obs on F2—10 for all AGN in the sample, (a— j) The individual data (crosses) are shown together with the 
binned data (filled circles, see Sec. 3.3 for details of binning). Errors on the binned data are plotted but are smaller than the size of the 
symbols, (k) Binned data of NGC 3516 resulting from the simple power-law model (filled circles) and power-law model modyfied by a 
Gaussian line and absorption edge with variable line and edge energies(open circles; models Ml and M6 in Fig. 4, see Sec. 3.3 for details). 
(1) Binned V — F2— 10 data of all AGN. Fluxes are normalized to the respective F2— 10, and we have subtracted the respective r o b s from 
the spectral slope values. NGC 5548 (filled squares) follows a different "slope vs. flux" relation than the other AGN (open circles). 
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source. For simplicity we will refere to it as the accretion 
rate hereafter. The errors on r o b s plotted in Fig. [5] are the 
standard deviation of r o t, s and we use them, instead of the 
standard error of the mean, to account for the scatter of the 
individual r o b s values about r o b s - 

The solid lines plotted in Figs. [S^i and [Sp indicate 
the average spectral slope of all AGN in our sample. In- 
dividual points appear to scatter randomly about this line. 
These plots suggest that the average spectral shape does 
not correlate with either the average luminosity (r = 0.2, 
Pnuli = 0.42) or BH mass (r = -0.38, P nuU = 0.13). On 
the other hand, panel (c) in Fig. [S] implies a positive cor- 
relation between F b s and mx.E (r = 0.51, P nu ii = 0.04). 
Sources with a higher accretion rate appear to have steeper 
average spectra as well. We found that a power-law model 
of the form y = ax 13 can fit well the data plotted in Fig. [5ji 
(X 2 = 8.14/8 d.o.f., Pnuii = 0.42). The best fitting parame- 
ters are: a = 2.7 ±0.3, and (3 = 0.08 ±0.02 (the error on the 
slope indicates that the correlation between r o b s and rhx,E 
is significant at the 4a level). 

3.2 Variations of the EW of the iron line 

An iron K a line is detected in many spectra of all sources 
(except for Akn 564, Fairall 9 and Mrk 766). To investigate 
the line's EW variations, we considered the spectral fit re- 
sults in the case when the line is detected at the 2a level in a 
spectrum of a source. The number of spectra which do show 
such a strong line, and the average equivalent width, EW, 
using the individual EW values from these spectra only, are 
listed in Tab. 0] The fourth column in Tab. H lists the \ 2 
values when we fit a constant to the EW light curves. The 
result is that we do not detect significant EW variations, 
except for NGC 5548, MCG -6-30-15, and NGC 3516. This 
is mainly due to the fact that, although the line is detected, 
the small exposure time of the individual RXTE observa- 
tions did not allow the determination of the line's EW with 
a high accuracy. We will present a more detailed study of 
the variability of the EW of the iron line in a forthcomming 
paper. 

Panels (a) , (b) and (c) in Fig. [6] show EW plotted as 
a function of L2-10, mx,E, and r o b s - Fig. |6p and Fig. [6ji 
imply that EW does not correlate strongly with either ac- 
cretion rate (r = -0.05, P nu ii = 0.88) or r obs (t = -0.14, 
PnuU = 0.65). On the other hand, Fig. [(Ji suggests that 
the line's EW anti-correlates with the source's luminos- 
ity: the strength of the line decreases with increasing X- 
ray luminosity (r = —0.52, P nu u = 0.09). This result is 
not statistically significant, most probably due to the small 
number of objects in the sample. However, this "equivalent 
width - luminosity" anti-correla tion that we observe (the so - 
called Iwasawa-Taniguchi effect. Ilwasawa fc TaniguchHll993l) 
is we ll established for nearby AGN (see e.g. iBianchi et al.l 
2007, and references therein). The solid line in Fig. [S^ indi- 
cates the best power-law fit to the data. The best-fit slope 
of -0.23 ± 0.06 is consi s tent w ith the slope of -0.17 ± 0.03 
found bv lBianchi et alj (|2007f ). 

3.3 The spectral variability of each source 

Crosses in panels (a-j) of Fig. [7] indicate the individual 
(r o b s , P2-10) points for each AGN in the sample. To bet- 



Table 4. The mean EW of the iron line and its standard devia- 
tion. N2 CT is the number of spectra with a 2<j detection of an iron 
line. 
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n 2ct 


EW 
(eV) 


x 2 


■fvar 

(%) 
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13 


682±351 


41 
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45 
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Figure 8. The binned "spectral slope - flux" data, as in Fig. [7] 
plotted with the best-fit function to the data plotted in Fig. [5}:. 
The sources were split arbitrarly between panels (a) and (b) for 
clarity of presentation. 

ter illustrate the long-term, "observed spectral slope - flux" 
relation in these objects, we grouped the (F b s ,P2-io) data 
in flux bins and calculated the (unweighted) mean flux and 
photon index in each bin. Each bin contains at least 20 mea- 
surements, and the number of bins is roughly equal for all 
sources. Filled circles in Fig. [7] indicate the (average spectral 
slope, average P2-10) points. Both the binned and un-binned 
data points in Fig. [7] suggest that, in almost all objects, the 
spectrum softens with increasing 2-10 keV flux. 

In Fig. [Tit we plot the binned "spectral slope - flux" 
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data for NGC 3516 using the best fit results of model M6 
(filled circles) together with those when we used the best 
fit results of model Ml (open circles; this is the model that 
shows the largest discrepancy in F b s with the statistically 
acceptable model M6 for this source). We fitted a power-law 
model to the Ml and M6 data plotted in this panel, and the 
best fit slopes differ by 0.08±0.02, with the Ml slope being 
flatter. Therefore, the exact shape of the T — -F2-10" plots 
in Fig. [7] does depend on the choice of the model used to fit 
the individual spectra of the sources. However, the overall 
trend stays the same in both cases: the spectrum steepens 
with increasing flux. 

In fact, the spectrum steepens with increasing flux in 
almost the same way in all sources. The opened circles 
in Fig. [7J indicate the binned T b s — .F2-10" data for all 
sources. We did not use different symbols for each source in 
order to emphasize the fact that, in all sources, the spec- 
tral shape varies with flux in a similar way. The only ex- 
ception is NGC 5548. Filled squares in Fig. indicate the 
'Tobs — i*2-io" data for this source. The spectral slope ap- 
pears to stay roughly constant although the 2-10 keV flux 
varies by almost an order of magnitude. A positive corre- 
lation between spectral slope and X-ray flux may exist at 
the highest luminosity bins, but there is also a hint of an 
anti-correlation at the lowest luminosity bins. 

The positive correlation between F b s and -F2-10 can be 
translated into a correlation between F b s and rhx,E, sim- 
ilarly to that between their average spectral slope and av- 
erage accretion rate (compare Fig. [5J;), but not necessarily 
the same. Figure[8]shows the binned (r o b s , Wx,e) data for all 
sources plotted together with the best-fit power-law model 
to the (r o b a , ?tix,e) data from Fig. [5}j. It is clear that, in at 
least some cases, such power-law cannot fit well the individ- 
ual "observed spectral slope - accretion rate" relations. 

To quantify the comparison between the average and 
the individual "observed spectral slope - accretion rate" re- 
lations, we fitted the (r o b s , t^x.e) data plotted in Fig. [8] 
with the same power-law model we used to fit the data 
shown in Fig. The model fits well the data of Akn 564 
(/?Akns64 = 0.10 ± 0.01, x 2 = 15/8 d.o.f., Pnull = 0.06), 
NGC 5506 (/W5506 = 0.12 ± 0.01, X 2 = 8.1/8 d.o.f., 
Pnuii = 0.43), Mrk 766 (/3 M rk7 66 = 0.16 ± 0.01, X 2 = 
6/8 d.o.f., p nul i = 0.75), and NGC 3516 (/Wssie = 
0.32 ± 0.02, x 2 = 5.8/7 d.o.f., p nu n = 0.56). Interestingly, 
/?Akn564 and /3ngc5506 are similar to the best fit slope in the 
case of the (F b s , »tix,e) plot. 



4 DISCUSSION 

The main aim of this work was to determine the shape of 
the X-ray continuum, and study its variability, in 10 X— 
ray luminous AGN which have been observed regularly with 
RXTE since 1996. We used a large number of spectra, which 
were taken regularly over a period of many years. As we 
argued in the introduction, since this time period is much 
larger than the characteristic time scale in these objects, it is 
possible that we have sampled the full range of their spectral 
variations. We therefore believe that we have determined, as 
accurately as it is possible at present, the "observed spectral 
slope - flux" relation for the AGN in the sample. In Section 
4.1 we discuss the implications from the average "observed 



spectral slope - accretion rate" correlation reported in Sec- 
tion 3.1. 

In Section 2.2.5 we argued that despite the limited spec- 
tral resolution of the PCA on board RXTE, the r o b s values 
we derived from the model fitting of the individual spec- 
tra of each source are representative of their actual X-ray 
spectral shape. However, this conclusion does not necessarily 
imply that r o b s = Fi ntr (Fintr being the intrinsic slope of the 
X-ray continuum). Due to the short exposure of the indi- 
vidual RXTE observations and the limited resolution of the 
PCA on board RXTE, we cannot constrain the properties 
of any warm absorbing material that may affect the spec- 
tra above 3 keV. Consequently, if there exists such a warm 
absorbing material in the vicinity of our X-ray sources, we 
would expect F b s < Tintr. Similarly, the possible presence 
of a significant reflection component (which we did not take 
into account in our models) will also affect the observed X- 
ray spectral shape. As a result, we would expect again the 
spectral shape of the continuum to appear harder than it 
actually is, i.e. r o b s < Fintr- 

Although all these mechanisms may operate in AGN, 
one of our main results is that the 'T b s — flux" relation is 
similar to all the objects in the sample, despite the fact that 
their BH mass and accretion rate are quite different. This 
common behaviour argues that their spectral variability is 
mainly driven by the same mechanism in all of them. In 
the following sections we demonstrate that, using various 
reasonable assumptions, the plots shown in Figs. [?H8]can be 
used to test various models for the AGN spectral variability 
that have been proposed in the last few years. 

4.1 The average X— ray spectral shape in AGN 
and Comptonization models 

One of the main results of our study is that the AGN in the 
sample do not have the same average spectral slope. Instead 
we found that F b s correlates positively with mx,E, that is 
AGN with higher accretion rate show softer X-ray spectra 
as well. This is in agreement with recent studies which also 
suggested a positive correla tion between spectral slope and 
accretion rate in AG N (e.g. IPorauet et aLll2004l ; lBianll2005l ; 
Sh cmmcr et all l2006t ISaez et all |200ST ) and Galactic black 



hole binaries ilWu fc Gull2008l . GBHs). 

iPapadakis et all (|2009h argued that the difference in the 
r o b s of the objects reflects a real difference in their intrinsic 
spectral slopes and it is not due to either the presence of a 
reflection component of different amplitude in the spectra 
of the sources or the effects of strong absorption in the 2-10 
keV band. Furthermore, the reality of the average "observed 
spectral slope - accretion rate" correlation we observed can 
also be justified from its similarity with the "F— accretion 
rate" relation that has been observe d in luminous quasar s 
(i.e. PG quasar luminosity and higher Ishemmer et al.ll2006h . 
since the energy spectrum in such sources is not expected 
to be significantly affected by either absorption or reflec- 
tion effects. We therefore believe that this correlation most 
probably reflects a "true/intrinsic" correlation between the 
intrinsic spectral slope and accretion rate in AGN. 

It is widely believed that hard X-rays from AGN are 
produced by Comptonization, i.e. by multiple up-scattering 
of seed soft photons by hot electrons in a corona located close 
to the black hole. The Comptonization process generally 
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produces power-law X-ray spectra. The main factor which 
determines the resulting X-ray spectrum is the so-called 
Compton amplification factor, A = (Ldiss + L a )/L a , where 
Ldiss is the power dissipated in the c orona and L B is the i nter- 
cepted soft luminosity. According t o Beloborodovl l|l999l ) and 
iMalzac. Beloborodov fc Poutanenl (120011 ) . F sa g(A - l)~ b , 
where o = 2.15-2.30, and b = 0.07-0.10 in the case of AGN, 
in which the energy of the input soft photons is of the order 
of a few eVs. 

If hard X-rays from AGN are indeed generated by 
the Comptonization process, then a possible explanation 
for the "spectral slope - accretion rate" correlation shown 
in Fig. [5h is that the ratio Ldi aa /L a in these objects de- 
creases proportionally with increasing accretion rate. In- 
deed, if Ldi aa /L a oc m _1 , then the r oc (A — l) _b relation 
implies that F oc rh b , with b = 0.07-0.1, as observed. 

The ratio Ldi aa /L a depends mostly on the geometry of 
the accretion flow, and the geometry may depend on the 
accretion rate. One possibility is that the cold dis c is dis- 
rupted in the inner region (e.g. lEsin et al.l Il997f ). As m 
increases, the inner radius of the disc may approach the 
radius of the innermost stable circular orbit, and so more 
soft photons are supplied to the hot inner flow causing 
an increase in L a . As a consequence, Ldi aa /L a should de- 
creas e and the X-ray spectrum should soften accordingly 
(e.g. IZdziarski. Lubiriski fc Smith! 1999V However, in this 
case, we will have to accept that the inner radius of the 
disc in AGN has not reached the innermost radius of the 
last circular stable orbit, even in the case of objects such as 
Ark 564 which probably accretes at almost its Eddington 
limit. 

Another possibility is that the coronal pla sma is moving 
away from the disc and emits beamed X-rays l|Malzac et al.l 
l200ll ). In this case, the Compton amplification factor de- 
pends on the plasma velocity, f3. As this velocity increases, A 
should also increase and the spectrum should harden. If this 
picture is correct, then our results imply that the plasma ve- 
locity should decrease with increasing accretion rate, which 
is somewhat opposite to what one would expect. Further- 
more, if the X-rays are beamed away from the disc, one 
would expect the EW of the iron line to decrease with in- 
creasing plasma velocity as well. Hence, one would expect 
the EW to increase with increasing accretion rate. However, 
we do not observe a significant correlation between the line's 
EW and accretion rate for the objects in the sample. 

The "observed spectral slope - accretion rate" relation 
could also be explained, quantitatively, if we assume that 
either the amount of reflection decreases or the effects of 
absorption become weaker as the accretion rate increases. 
In the first case, we should also expect to detect an anti- 
correlation between EW and accretion rate, which is not 
the case. In the second case, if the warm absorber is located 
far away from the central source and the weakening of the 
absorption is due to the decrease in its opacity, by increasing 
its ionization parameter for example, we should expect r o b s 
to correlate with the source luminosity, which is not the 
case. It is more difficult to anticipate the "observed spectral 
slope - accretion rate" relation in the case when the absorber 
originates from an outflowing disc wind. The reason is that 
no physical models predict at the moment how does the 
accretion rate control the mass loss rate and the ionization 
state of the wind, and hence its opacity. 
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Figure 9. The 'T t, s — F2— 10" relation for three r; ntr values in 
the case the X— ray continuum varies only in flux, and there exists 
absorbing material whose ionization parameter varies proportion- 
ally with the continuum flux ( solid and dashed lines correspond 
to N H = 5 X 10 22 and N H = 10 23 cm -2 , respectively). 



4.2 The spectral variations within each source 

4-2.1 Intrinsic spectral slope variability 

One possible explanation for the "observed spectral slope 
- flux" relations shown in Fig. [7] is that the observed flux 
variations correspond to m variations in each source. We 
argued above that it is possible for the Compton amplifica- 
tion factor, A, to decrease proportionally with increasing m 
for different AGNs. Similarly, this factor may vary with rh 
also in individual objects, Iin which case we would expect 
r oc -F|_ 10 , where b ~ 0.07-0.1 (see discussion in the pre- 
vious section). However, we found that a power-law model 
fits the 'T— F2-10" relation of just 4 AGN in the sample. 
Furthermore, only in the case of Ark 564 and NGC 5506 
the best-fitting slope is consistent with the expected value 
of ~ 0.07-0.1. Therefore, the flux related spectral variations, 
at least in the other objects, most probably have a different 
origin. 



4-2.2 The case of variable absorption. 

Spectral variations can also be caused by variations in 
the column density, covering fraction (say /), and/or the 
opacity of an absorber while Tjntr r emains constant (e.g. 
iTurner et afll2007t iMiller et al1l2008l ). It is rather unlikely 
that variations of either the column density or the covering 
factor are the main driver for the "spectral slope - flux" re- 
lation. Under this hypothesis, the Nh or / variations flatten 
the observed spectral slope and and also reduce the con- 
tinuum flux, as observed. In other words, in this case the 
observed flux variations in AGN should not be intrinsic, but 
they should be caused by the variations in the column den- 
sity or the covering factor of the absorber. However, many 
sources in the sample show significant flux variations on time 
scales as short as a few hundred seconds. Therefore, the con- 
tinuum flux variability in these objects must be intrinsic to 
a large extent, as neither Nh nor / are expected to vary on 
these time scales. 

On the other hand, it is possible that the variability of 
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the ionization state of the absorber, due to intrinsic varia- 
tions of the continuum flux, can affect the spectral shape of 
the continuum and result in "spectral slope - flux" relations 
which are similar to the relations we observed: a decrease 
in the continuum flux will decrease the ionization state of 
the absorber, increasing its opacity and hence flattening the 
observed spectrum. 

To investigate this possibility we considered the 
absori*powerlaw model in XSPEC, where absori is an 
ionised absorption model (jPone et al.l ll992). The main pa- 
rameters of this model are the ionisation parameter £ and 
the column density, JVh, of the absorbing medium, as well as 
Tintr and Npl, normalisation of the power-law component. 
We created model spectra assuming Fi nt r = 1.5, 2 and 2.5, 
and various values of Npl (from 1 up to 5000). We then con- 
sidered three values of column density, Nh = 10 22 , 5 x 10 22 
and 10 23 cm -2 and, for each I\ n tr, we assumed that the ion- 
isation parameter of the absorber is changing proportionally 
with N PL (i.e. £ txN PL ). 

We then used XSPEC to fit a simple PL model to the re- 
sulting spectra and determined r o b s and F2-io- In Fig.[9]we 
plot r o b s as a function of the observed flux for the three Tintr 
values we considered and in the case when Nh = 5 x 10 22 
and Nh = 10 23 cm -2 (solid and dashed lines, respectively; 
for simplicity, we assumed, arbitrarily, that £ =Npl). As 
expected, the observed spectral slope flattens as the source 
luminosity decreases. The r o b a change is more pronounced 
in the case of the absorber with the larger absorbing column 
(the r o b s variations are of low amplitude in the case of the 
Nh = 10 22 absorber, and for this reason we do not plot the 
respective models in Fig. [9j . 

We fitted the curves shown in Fig. [9] to the data plot- 
ted in Fig. [7] by allowing them to shift both in the vertical 
(i.e. j/ m odci = r o bs + 2/ s hift), an d the horizontal direction (i.e. 
^modci = F2-10 x s s hift). These model curves fitted well the 
"F obs - flux" data of NGC 5548 only. A N H = 5 x 10 
cm 

*~ 

the spectral variations we observed in this object. In all other 
cases, the observed "spectral slope - flux" relations are sig- 
nificantly steeper than any of the model curves plotted in 
Fig.H 

We conclude that opacity variations, due to changes in 
the continuum flux, of a single absorber, whose column den- 
sity remains constant, cannot explain the observed spectral 
variations in most AGN in the sample. Of course, if there are 
multiple layers of absorbing material, the resulting "spectral 
slope - flux" model curves will be different. However, it is not 
easy to model such a situation without some apriori knowl- 
edge of how many absorbing layers there may be, what their 
average ionization is, etc. 



4-2.3 The case of a constant reflection component 

It is possible that the observed spectral variations are due to 
the combination of a highly variable (in flux) power-law con- 
tinuum (with ri n tr=constant) and a constant reflection com- 
ponent (e.g. iTavlor Uttlev fc M c Hardvl 120031 ; iPonti et al.1 
l200d : lMiniutti et al.ll2007n ~ In this C£LS6, ELS the amplitude of 
the continuum increases, the relative contribution of the re- 
flection component to the observed spectrum will decrease, 
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2 absorber whose ionisation parameter varies between 
30 — 180, can explain satisfactorily (y 2 = 14.5/8 d.o.f.) 
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Figure 10. r o b s plotted as a function of -F2— 10 (normalized to the 
mean model flux) in the case of a power law with ri ntr = 2 which 
varies only in normalization and a constant reflection component 
in the case when the reflecting material is neutral (i.e. log(£) = 0; 
solid line) and ionized (log(£) = 3; dashed line). The dotted line 
indicates the "spectral slope — flux" relation in the case when ri ntr 
varies with flux (r; ntr oc -F^— °lo) an< ^ there also exists a constant 
reflection component (with R = 1). 



and r o bs will increase, reaching Tintr when the continuum 
amplitude is large enough. 

To investigate this possibility, we considered a 
"power-law plus constant reflection" model. To calculate 
the reflected component we used the pexriv code by 
iMagdziarz fc Zdziarskil(|l995l ). available in XSPEC. The main 
model parameters are: (i) the X-ray photon index of the in- 
trinsic power-law, Fi nt r, (ii) the ionisation parameter of the 
reflecting medium, £, and (iii) the amplitude of reflection, 
R. We computed model spectra assuming Fi ntr = 2, various 
values of the power law normalization, Npl (from 3 x 10 -5 
up to 0.1), while the normalization of the reflected compo- 
nent, N ren , was kept fixed (the pexriv model was defined, 
arbitrarily, so that \R\ = 0.1 and N rGn = 0.1, where iV rc fl 
means normalization of the power-law continuum that is re- 
flected). Since the power-law was defined explicitly in our 
model, we used pexriv with R < 0, which outputs only the 
reflected component. 

We then used XSPEC to fit a power-law to the model 
spectra in the 2-10 keV band and determined model r o b s - 
In Fig. [10] we plot Fobs as a function of the total 2-10 keV 
flux (normalized to the mean model flux) in the case when 
the reflecting material is neutral (i.e. log(£) = 0; solid line) 
and ionized (log(£) = 3; dashed line), and the slope of the 
power law is fixed at Ti ntI = 2. As expected, as the flux 
decreases, the relative strength of the reflecting component 
increases, and r o bs becomes harder. 

The shape of the model curves in the case of different 
Fintr values is identical to the ones shown in Fig. \TU\ ex- 
cept for the fact that, in the high flux limit, they saturate 
at their respective Fir.tr value. Similarly, the shape of the 
model curves remains the same if we increase the total flux 
while keeping the relative strength of the power-law and the 
reflection component (and hence F b s ) the same. As a result 
the curves will only shift to a higher flux. 

We fitted the curves shown in Fig. [10] to the data plot- 
ted in Fig. [7J by allowing them to shift both in the vertical 
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(i.e. j/model = r obs + 2/shift), and the horizontal direction (i.e. 

-F2-10, mod x Xahift; ^2-io, mod is the model 2-10 keV 
flux). The shift along the y— axis is necessary because Tintr 
may not be equal to 2 in all objects. The best fit vertical 
offset value, j/bf, will determine the best fit intrinsic slope 
for each object, r^tr = 2 + y hi - On the other hand, the hori- 
zontal offset shift is necessary, given the different flux levels 
of the objects, and the fact that the strength of the reflec- 
tion component in each object, with respect to that of the 
power- law continuum (i.e. the average R in each object, R), 
is not known apriori. In fact, the best fit horizontal offset 
value, Xbi, can be used to determine R as follows. 

The reflection amplitude, R, determines the fraction of 
the X-ray photons that are reflected by the disc, so that 
R oc -F r es /-Fpl (where -Fred and -Fpl are the fluxes of the re- 
flected component and continuum, respectively). In our case, 
the average source flux, i ? 2-io,obs, corresponds to a model 
PL flux of FpL.mod = (-F^-io.obs/^bf) --F r cfl,mod- The model 
curve was constructed by assuming, arbitrarily, the constant 
reflection component when R = 0.1 and Npl = 0.1, there- 
fore R = 0.1 x F PL , mod (0.1)/F PLimod , where F PL ,mod(0.1) is 
the 2-10 keV flux of the power law with Npl = 0.1. 

The model curves plotted in Fig. [I0]fitted well the data 
of all objects, except for NGC 3227 {\ 2 = 55.1/8 d.o.f.), 
NGC 5506 (x2 = 22.2/8 d.o.f.), and NGC 5548 (x2 = 
40.4/8 d.o.f.). In the case of MCG -6-30-15 and NGC 3516 
it was reflection from ionized material that fitted the data 
best. We plot the 'T bs ~ F2_io" data for the objects in the 
sample in Fig. 1111 Best fit lines are indicated by the solid 
lines. 

The best model fits do not provide a statistically accept- 
able fit to the data of NGC 3227, NGC 5506, and NGC 5548. 
In the case of NGC 5548, the model describes rather well the 
"Fobs— flux" data, except for the lowest flux point (the best 
fit model line plotted in Fig. [TT] for NGC 5548 corresponds 
to reflection from ionized material). The disagreement be- 
tween best model fit and data points is more pronounced 
in the case of NGC 3227, especially at low fluxes. This is a 
source where a possible absorption event, which lasted for 
a few months, has been observed in its RXTE light curve 
(|Lamer. Uttlev fc M c Hardvll2003l 1. In this case, one would 
expect r o b s to be harder than the model predictions at a 
given (low) flux, which is opposite to what is observed. How- 
ever, such an event should also shift the model line to lower 
fluxes. Therefore it is difficult to reach final conclusions re- 
garding the quality of the model fit to the NGC 3227 data. 

The best fit intrinsic spectral slopes ranged from Tmtr ~ 
1.8 in NGC 5548, to r intr ~ 2.8 in Ark 564. Furthermore, 
the best fit Sbf values imply large average R values: from 
R ~ 0.5 in the case of NGC 5548 up to R ~ 15 in the case 
of NGC 3516. The median R is 4.7. We did not detect any 
correlation between the resulting best-fit Fi ntr and R values. 

The reflection amplitude R is essentially equal to the 
solid angle covered by the reflecting material and is equal 
to 1 in the case of a point source on top of a disc at large 
heights (i.e. in the case when Q = 2n). Values as large as 
R ~ 5 suggest either a peculiar geometry (which is diffi- 
cult to envisage) or non-isotropic emission from the X-ray 
continuum in which case the X-ray flux on the reflecting 
material is muc h larger than the flux wh ich reaches the ob- 
server (see e.g. [Miniutti fc Fa bian 2004j, for the case of an 
X-ray source which is located very close to a rapidly rotat- 
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Figure 11. The r o b s — F2-10 relations for the objects in the 
sample together with the best fit curves for the models described 
in sections 4.2.3 (solid line) and 4.2.4 (dashed line). 



ing black hole, in which case relativistic light bending effects 
are important). 



4-2.4 The case ofT ln t r variations plus a constant 
reflection component 

We considered next the possibility that the X-ray con- 
tinuum in AGN consists of a constant reflection component 
and a power-law continuum with variable Npl and Tintr- 
This is a combination of the two cases we considered above: 
X-rays are produced by thermal Comptonization, the flux 
variations (i.e. the Npl variations) reflect intrinsic varia- 
tions in the accretion rate, and Tintr = 2.7m ' 08 (this is the 
best power-law model fit to the Tintr — rh relation plotted in 
Fig. [SJi) . The X-ray source illuminates a reflecting material, 
which effectively responds to the average X-ray continuum 
(i.e. to a power-law of ~ Fi ntr and Npl). This could be for 
example the case of a reflector which is located away from 
the X-ray source. 

To investigate this possibility, we considered again the 
model pexriv+powerlaw in XSPEC. This time, we created 
model spectra assuming 8 values of the power law nor- 
malization (from 10~ up to 0.1) and a spectral slope of 
Fintr = 2.7N p 'l 8 (the arbitrary choice m =Npl was made 
for simplicity; it does not affect our results, as long as the 
resulting model curves are allowed to shift freely along the 
x — axis when we fit them to the data). Then, for each pair 
of values (Npl, Tintr) we considered the respective reflection 
component (from neutral material) in the case when R = 1, 
and we added it to all 8 power-law spectra. In this way, we 
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constructed 8 sets of model curves (8 model spectra in each 
set). 

We then used XSPEC to fit a simple PL model to the 
8 spectra of each model set and determined F b s and the 
spectrum flux in the 2-fO keV band (i^-io)- The dotted 
line in Fig. \TU\ indicates r obs in the case when we added to 
all 8 spectra the reflection component (with R = 1) from 
a Tintr = 2 power-law continuum, plotted as a function of 
F2-10 (normalized to the mean model flux). 

We fitted the 8 T b a - -F2-10" model curves to the 
observed "spectral slope - luminosity" relations plotted in 
Fig. [7] by allowing them to shift both in the vertical (i.e. 
2/modei = r o bs + J/shift), and the horizontal direction (i.e. 
a^modoi = x S Mtt x -F2-io,mod). These model curves fitted well 
the 'T — rh" data of all objects, except for NGC 4051 and 
NGC 5548. 

The dashed lines in Fig. [11] indicate the best fit 
curves for this model. Interestingly, the model fits well the 
NGC 3227 data, a result which implies that the flat r o t, a 
at low fluxes may correspond to intrinsically flatter spectra, 
and not an increase in the external absorption by cold gas. In 
the case of NGC 4051, the fit is not statistically acceptable 
(X 2 = 20.9/8 d.o.f., pnuii = 7.4 x 10~ 3 ) mainly because of 
the lowest flux point. Furthermore, the model fails entirely 
to describe the NGC 5548 data. This result is obvious given 
the fact that the model curves in this case have no flat part 
that could fit, even approximately, the observed NGC 5548 
data. 

In principle, we should not allow a y— axis shift of the 
model curves in this case. The intrinsic spectral slope is 
uniquely determined by the Fi n t r — Npl relation we assumed. 
However, the assumed relation was determined by the PL 
model fit to the "average spectral slope - average accre- 
tion rate" data plotted in Fig. [SJj, hence there is an un- 
certainty associated with its normalization and slope. The 
average best fit y— axis shift for all objects is 0.13±0.08, i.e. 
it is consistent with zero. This result implies that the rela- 
tion Fintr ~ 2.7m 08 is in agreement with both the average 
"spectral slope - accretion rate" relation of all objects and 
with the individual "slope - luminosity" relations of each 
AGN. 

The Ark 564 data are fitted best in the case when the 
constant reflection component is the one which corresponds 
to r in tr = 2.45. The Mrk 766 and NGC 4051 data are_fitted 
best in the case of constant reflection component for r in tr = 
2.25 and 1.9, respectively. In all other cases, the addition of 
a constant component for Tintr = 1.7 is needed to fit the 
data best. 



5 CONCLUSIONS 

We analyzed more than 7,500 RXTE PCA data of 10 nearby 
AGN with the goal of studying their spectral variability 
properties. The objects were observed regularly with RXTE 
over many years since 1996. Given the large number of obser- 
vations, spread over many years, we were able to determine 
accurately their average observed spectral slope, luminos- 
ity and accretion rate. Using the best model-fit results from 
the hundreds of spectra of each source, we also determined 
their long-term "observed spectral slope - X-ray flux" re- 
lation. This kind of long term variability analysis can not 



be accomplished with the rather short AGN observations 
currently performed by Chandra, XMM-Newton or Suzaku, 
despite their higher effective area and superior spectral res- 
olution. Our main results are as follows: 



1. The average spectral slope is not the same in all objects. 
It does not correlate either with BH mass or luminosity. 
However, Y correlates positively with the mass accretion 
rate: objects with higher accretion rate have steeper spec- 
tra, and T ~ 2.7m 08 . This result can be explained if the 
Compton amplification factor, A, decreases proportionally 
with the accretion rate in AGN. 

2. We detected the iron Ka line in many spectra of seven 
sources in the sample. The average EW of the line anti- 
correlates with the average 2-10 keV luminosity (i.e. the 
so-called Iwasawa-Taniguchi effect). Our results are in agree- 
ment with recent results which are based on recent studies 
of large samples of AGN l|Bianchi et al.ll2007T ). 

3. The spectra of each source become steeper with increas- 
ing flux. This is a well known result. What we have shown 
in this work that this trend holds over many years, and over 
almost the full range of flux variations that AGN exhibit. 
We also found that the r oDS — F2-10 relations are similar 
in all objects, except for NGC 5548. There are many rea- 
sons which can cause apparent and/or intrinsic variations 
in the X-ray spectrum of AGN. It is possible that they all 
operate, to a different extent in various objects. However, 
the common spectral variability pattern we detected implies 
that just a single mechanism is responsible for the bulk of 
the observed spectral variations. 

4. NGC 5548 displays limited spectral variations for its 
flux variability. Although uncommon, this behaviour has 
also been observed in other Typ e I Seyfert, PG 0804+761 
(jPapadakis. Reig fc Nandrall2003l) . This behaviour, different 
to what is observed in most (Type I) AGN, raises the issue 
of different spectral states in AGN, just like in GBHs. In- 
vestigation of this issue is beyond the scope of the present 
work. 

5. The scenario in which the spectral variability is caused by 
absorption of X-rays by a single medium whose ionization 
parameter varies proportionally to the continuum flux vari- 
ations fails to account for the observed spectral variations. 

6. A "power-law, with constant Ti n tr and variable flux, 
plus a constant reflection component" can fit the observed 
'T — F2-10" relations of most (but not all) objects in the 
sample. The flux of the reflection component necessary to 
explain the observed spectral variability is rather large, im- 
plying average reflection amplitudes of the order of R ~ 5. If 
true, this result implies that the primary source of X-rays is 
located close to a maximally rotating Kerr black hole, where 
relativistic effects (such as light bending) are expected to be 
strong. 

7. The observed 'T — J2-10" relations (except for 
NGC 5548) can be explained if we assume that the power- 
law continuum varies both in flux and shape as Fintr oc i^-To 
(as implied by the average "spectral slope - accretion rate" 
relation that we observed) and the reflecting material re- 
sponds to the average continuum spectrum, with R = 1. 
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